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P-i\ ABSTRACT 

Q . GRB 061007 was the brightest gamma-ray burst (GRB) to be detected by Swift and 

$H ' was accompanied by an exceptionaUy luminous afterglow that had a y-band magni- 

I tude < 11.1 at 80 s after the prompt emission. From the start of the Swift observations 

d • the afterglow decayed as a power law with a slope of ax = 1.66±0.01 in the X-ray and 

I aopt = 1-64 ± 0.01 in the UV/optical, up to the point that it was no longer detected 

• i-H , above background in the optical or X-ray bands. The brightness of this GRB and the 

' similarity in the decay rate of the X-ray, optical and 7-ray emission from 100 s after the 

H I trigger distinguish this burst from others and present a challenge to the fireball model. 

. 5t 1 The lack of a cooling or jet break in the afterglow up to ^ 10^ s constrains any model 

that can produce the large luminosity observed in GRB 061007, which we found to 
require either an excessively large kinetic energy or highly collimated outflow. Analysis 
of the multi-wavelength spectral and high-resolution temporal data taken with Swift 
suggest an early time jet-break to be a more plausible scenario than a highly energetic 
GRB. This must have occurred within 80 s of the prompt emission, which places an 
upper limit on the jet opening angle of 9j = 0.8°. Such a highly collimated outflow 
resolves the energy budget problem presented in a spherical emission model, reducing 
the isotropic-equivalent energy of this burst to E!y°^^ — 10^" erg, consistent with other 
GRBs. 
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1 INTRODUCTION 

During their initial outburst, gamma-ray bursts (GRBs) are 
the most luminous source of electromagnetic radiation in the 
Universe, releasing isotropic-equivalent energies of the order 
of 10''*-10^* erg in 7-rays alone on timescales of 10~^-10^ s. 
The favoured mechanism to produce such a large amount 
of energy on such a short timescale involves accretion onto 
a newly formed black hole produced either in a merger of 



two compact objects or from the gravitational core collapse 
of a massive star. There appear to be two distinct popula- 
tions of GRBs, distinguished by the durat ion of their prompt 
emiss ion and their spectral properties (|Kouveliotou et al] 
Il993h . whereby long GRBs (prompt emissions lasting > 
2 s ^) are typically softer than short GRBs (< 2 s) 
(|Sakamoto et al.l |2006| ). The mos t robust mode l for long, 
soft GRBs is the coUapsar model (|Wooslevlll993l ). in which 
a GRB results from the energy released during the death 
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of a massive star. Evidence supporting this model includes 
the large number of lon g GRBs located in regions of ac- 
tive star-formation (e.g. iTanvir et al.l l2004l : iFruchter et al.l 
I2OO6I ). and the association of several near by GRBs with 
super novae, e.g. GRB 980425/S N1998bw jKulkarni et al] 
1 19981 ). GRB 030329 /S N2003dh jHiorth et all |2003| ) and 
GRB 060218/SN2006aj (|Campana et al.ll2006l ) 

There is evidence to suggest that collimated outflows re- 
duce the amount of energy emitted by a GRB from isotropic 
emission by ~ 2 orders of magnitude, where the beam- 
corrected energy is clustered at lO^^-lO^^ erg (Frail et al. 
2001; Bloom, Frail & Kulkarni 2003). Bipolar jets are a nat- 
ural consequence of accretion onto a black hole and are ob- 
served in many high-energy astrophysical phenomena (e.g. 
AGN, pulsars, micro-quasars), where the ejecta take the eas- 
iest exit path, along the rotational axis. 

A signature of highly relativistic jets is an achromatic 
break in the broadband afterglow light curves, caused by 
the expansion of the radiation beaming angle as the bulk 
Lorentz factor, F, decreases (Rhoads 1999; Sari, Piran & 
Halpern 1999). Once the Lorentz factor has decreased such 
that 6j = l/F, where 6j is the half jet opening angle, the rate 
of decay increases, resulting in a break in the light curve. 

One of the scientific aims of the GRB dedicated mis- 
sion, Swift l|Gehrels et al.ll2004l ). was to pin down the range 
of coUimation present in GRBs. Swift is a rapid-response 
mission with multi-wavelength capabilities designed specif- 
ically to study the early time evolution of GRBs and their 
afterglows. The spacecraft is equi pped with three telesco pes; 
the Burst Alert Telescope (BAT; iBarthelmv et allbOOSl ) has 
a 15-150 keV energy range with a partially coded response 
up to 350 keV, the X-Ray Telescope (XRT; iBurrows et al.l 
I2OO5I ) which covers the 0.3-10 keV energy r ange, and the 
Ultra - Violet and Optical Telescope (UVOT; iRoming et al.l 
I2OO5I ). which can observe in six colours and a broadband 
white band filter, covering the wavelength range from 1600 A 
to 6500 A. 

Swiffs slewing capabilities provide GRB afterglow light 
curves from ~ 60 s after the trigger, and well-sampled X-ray 
and optical light curves have revealed far fewer achromatic 
breaks than had been expected prior to launch (Ghirlanda, 
Ghisellini & Lazzati 2004; O'Brian et al. 2006). This may 
be an indication that jet breaks occur later than previously 
believed, at a time when the afterglow is no longer detected 
above the background level. In the X-ray band this is at 
around 10® s after the prompt emission. The time of the 
jet break, tj, increases with the jet opening angle, 6j, as 
tj oc (n/ Ek)~^^^ erg, for a constant density circumburst 
medium with density n and GRB kinetic energy Ek- Large 
average jet opening angles, small circumburst densities or 
large GRB kinetic energies may, therefore, all contribute to 
delaying the time of the jet-break to when the afterglow is 
no longer detectable by Swift, although tj is only weakly 
dependent on the latter two parameters listed. 



^ The duration of the prompt emission for the BATSE GRB 
sample showed a bimoda l distribution with separation at ~ 2 s 
llKouveliotou et al.|[l993l') . However, in the light of Swift the dis- 
tinction betw een short and long GRBs has become les s clear (e.g. 
GRB 060614: [Gehrels et al.ir2006l : IZhang et al.ll2007d) . and both 
temporal and spectral properties need to be considered to distin- 
guish between these two populations of GRBs. 



GRB 061007 is an example of a burst that showed no 
evidence of a break in the X-ray or UV/optical afterglow 
from 80 s after the prompt emission, until it was no longer 
detected above background, at ~ 10® s. Prompt observa- 
tions by ground based observatories measured the redshift 
of this burst to be z=1.26 (Osip, Chen & Prochaska 2006; 
Jakobsson et al. 2006). At this redshift, GRB 061007 had 
one of the highest isotropic-equivalent energies ever seen, 
releasing ~ lO^** erg in 7-radiation alone, wi th an optical 
afterglow comparable to that of GRB 990123 (|Akerlof et al.l 
Il999h . with V < 11.1 mag at 71.7 s after the BAT trigger. 

In this paper we describe in detail the 7-ray, X-ray and 
optical observations taken by Swift. From the high resolu- 
tion early time data we constrain the type of environment 
in which the burst occurred and determine the conditions 
that distinguish this burst from others. In section [2] we de- 
scribe the data reduction and analysis techniques used, fol- 
lowed by an analysis on the possible physical parameters and 
emission mechanisms that could explain the observations in 
section 13.1113.31 In these subsections we investigate both a 
wind and ISM circumburst environment, as well as spherical- 
equivalent emission and collimated emission. We also use the 
afterglow spectral energy distribution (SED) to analysis the 
properties of the circumburst medium of GRB 061007 in 
section 13.41 Our results are discussed in section 13.51 and in 
section section|4]we summarise our conclusions. Throughout 
the paper all errors quoted are la unless otherwise specified 
and the temporal decay rates, a, and spectral indices, /3, 
are defined such that Ft oc t~" and oc , respectively. 
Luminosities and the energy release are calculated assuming 
a standard cosmology model with Hq = 70 km s~^ Mpc~^, 
fiif = 0.3 and Q.x = 0.7. 



2 OBSERVATIONS AND ANALYSIS 
2.1 BAT 

The BAT on board Swift triggered on GRB 061007 on 7*'' 
October, 2006 at 10:08:08 UT, resulting in an immediate 
slew to point the narrow field instruments in the direction 
of the burst. From here on we shall refer to the time of 
the trigger as T. The refined BAT position for GRB 061007 
is RA = 03''05"'11.8', Dec = -50°29'47.4" (J2000) with 
an uncertaint y (systematic and stat istical) of 0.9' at 90% 
containment (iMarkwardt et al]|2006l ). 

The BAT light curve consists of three large and distinct 
peaks and is shown in Fig. [1] The time-interval over which 
90% of the 15-150 keV 7-radiation is emitted was T90 = 
75 ± 5 s, although emission above the background level was 
detected out to ~T-|-500 s, which is most pronounced in 
the soft energy band (t op panel of Fig. [TJ. The time lag 
(jNorris fc Bonnellll2006l ) between the peak emission in the 
50-100 keV and the 15-25 keV energy bands was 16lg ms, 
and between the 100-350 keV and 25-50 keV bands the lag 
was 21^1 ms. 



2.2 XRT 

The first XRT observation was a centroiding image taken 
80.4 s after the BAT trigger, and a bright uncata- 
logued source was found at RA = 03''05'"19.5% Dec = 
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Figure 1. Prompt emission light curve of GRB 061007 shown in 
the four BAT energy bands. 



break in tlie liglit curve, resulting from either the edge of 
the jet coming into the line of sight or the migration of the 
cooling frequency, Vc, through the observing energy bands. 
The cooling frequency defines the frequency above which 
electrons no longer lose a significant fraction of their energy 
to radiation. The cooling frequency is expected to evolve as 
t~^^'^ for an ISM-like circumburst environment (Sari, Piran 
& Narayan 1998), and as t^^'^ for a wind-like circumburst en- 
vironment (I Chevalier &: 111120001) ■ and as it migrates through 
the observing energy band the drop in radiation at Vobs > 
produces a break in the light curve. 

To determine the 3(t lower limit on the time at which a 
jet or cooling break could have occurred, we applied a broken 
power law fit modelled on the cooling break, for which 02 = 
a.\ -t 0.25 dSari et al.lll998l ). and on a jet break, where 02 = 
Qi -t 0.75 jRhoaddligQgi r and in each case determined the 
earliest time that the break could occur within the observing 
time interval without deviating the goodness of fit by more 
than 3a from the best-fit power law model. Our analysis 
indicates that a cooling break could not have occurred earlier 
than 10^ s after the BAT trigger to 99% confidence, and any 
jet break after the onset of the afterglow observations must 
have happened at tj > 2.65 x 10^ s. 

A power law fit to the early WT spectrum (<T-t600 s), 
where the Galactic column density is fixed at A'^^f = 2.13 x 
10 cm' (jPickev fc Lockman|[l990l ), gives a best-fit host 
galaxy column density of AfH((WT)) = (5.3 ± 0.2) x 10^^ 
cm^'^ at a redshift oi z — 1.26 for a spectral index 
/3x((WT)) = 0.96 ± 0.01. Spectral analysis of the PC data 
during the time interval T-l-5000 s to T-f 25,000 s gives a 
best fit host galaxy column density Nh (PC) = (5.1±0.6) x 
10^^ cm"^, and spectral index /3x((PC)) = 0.92±0.5, consis- 
tent with the WT data, and therefore indicating that there 
is no spectral evolution in the X-ray data. 



-50°30'01.9" wit h a 90% containment radius of 3.5" 
IVetere et al.ll2006l ). 

At 86.6 s after the BAT trigger the XRT began tak- 
ing data in Window Timing mode (WT), in which data are 
stored with a 1.7 ms time resolution and 1-dimensional imag- 
ing. At T-l-2000 s the count rate fell to below ~ 1 count s~^ 
and the automated sequence switched to Photon Counting 
(PC) mode, in which the time resolution is 2.5 ms and both 
imaging and spectroscopic information are available. The 
X-ray afterglow was detected for 10® s after the BAT trig- 
ger before it fell below the detection threshold. During the 
first 2200 s of observations the count rate of the burst was 
high enough to cause pile-up in both WT and PC mode 
data. To account for this effect, the WT data were extracted 
in a rectangular 40x20-pixel region with a 9 x 20-pixel re- 
gion excluded from its centre. The size of the exclusion re- 
gi on was determined following the procedure illustrated in 
[Romano et al.l l|2006ah . The standard correc tion for pile-up 
was ap plied to the PC data, as described bv lVauehan et aP 
()2006l ). Due to a large amount of pile- up in the first 30 s of 
WT observations, these data were not used in our spectral 
analysis. 

The 0.3-10 keV light curve is well fit by a power law 
with decay index a = 1.66 ± 0.01 from T-l-80 s until at 
least T-l-10® s, when it reached the detectability threshold 
of XRT (see Fig. [Sjl . We also considered the presence of a 



2.3 UVOT 

The UVOT began observing GRB 061007 71.7 s after the 
BAT trigger, at which point it took a 9 s settling exposure 
in the V-filter. This was followed by a 100 s exposure in the 
white light filter (A ~ 1600-6500 A) and a 400 s exposure in 
the l/-band. After this the automated sequence rotated twice 
through the UVOT filters taking a series of short exposures 
(y,VVM2,UVW\,U,B, white, UVW2) (10 s for B and 
white filters and 20 s for the rest). The central wavelengths 
of the optical and UV filters are 5460 A, 4370 A 3460 A, 
2600 A, 2260 A and 2030 A for the V, B, U, UVWl, UVM2 
and UVW2, respectively. A further 400 s V-band exposure 
and 100 s white exposure were taken, followed by a further 
series of rapid rotations around the filter wheel up to 2260 s 
after the trigger. With the exception of the second 100 s 
white band exposure, all observations up to this time were 
taken in event mode, which has a time resolution of 11 ms, 
and after this image mode exposures were taken. The re- 
fined position of the afterglow is RA = 03''05"'19.6^ Dec 
= -50°30'02.4" to a certainty of 0.5" (J2000). The Galac- 
tic extinction along this line of sight is E{B — V) = 0.021 
(Schlegel, Finkbeiner & Davis 1998). 

For all short exposures (10 s or 20 s) and exposures 
taken later than T-l-2260 s photometric measurements were 
obtained from the UVOT data using a circular source ex- 
traction region with a 6" radius for the V , B and U optical 
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Figure 2. GRB afterglow light curve shown in all six UVOT lenticular filters, in the X-ray (0.3-10 keV) and in the 7-ray (15—350 keV). 
All light curves have been shifted vertically with respect to the V band filter for clarity. The B, U, UVWl, UVM2, UVW2 and white 
filter have all been normalised to the V band light curve. The left hand axis corresponds to the V-band light curve, and the count rate 
axis on the right applies to the X-ray and 7-ray light curves. The best fit decay indices to the combined UVOT light curve, the X-ray 
light curve and the BAT light curve are ajjYOT = 1.64 ± 0.01, axRT = 1.66 ± 0.01 and qbat = 1-54 ± 0.12, respectively. Error bars 
are smaller than the symbol size if not visible. 
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filters and a 12" radius for tfie three UVOT ultra-violet fil- 
ters to remain compatible with the current effective area and 
zeropoint calibrations^. The background was measured in 
a source-free region near to the target using an extraction 
radius of 20". 

The first V-band settling exposure and first white-band 
exposure were affected by coincidence loss to a degree that 
is outside the photometrically calibrated range of UVOT; 
for a V^-band settling exposure complete saturation within a 
standard 6" aperture implies that this afterglow is brighter 
than V = 11.1 for the whole of the 9 s exposure. 

During the first white band exposure the source is so 
bright that relative photometry can be obtained using a 
very large aperture, which includes the wings of the point 
spread function, where coincidence loss is not a problem. 
The exposure was divided into 5 s bins, and source counts 
were obtained from an extraction region of radius ~ 26". To 
calibrate this photometry, observations of the white dwarf 
standards GD50 and HZ2 were taken with UVOT and the 
white filter on October 18"" 2006. From these data we deter- 
mined the relationship between observed count rate within 
the 26" aperture and corrected count rate within a standard 
6" UVOT source aperture, which we assume to be linear. To 
apply this calibration to GRB 061007 we apply constant off- 
set to this relation to account for the different background 
count rate between the GRB 061007 observations and the 
white dwarf standard count rate observations. This offset 
was chosen to match the mean white-band count rate in the 
saturated exposure to the white filter light curve extrapo- 
lated from later times; we caution that calibration of this 
first white exposure is crude and thus the apparent steepen- 
ing of the light curve towards the end of the exposure may 
not be real. However, the overall behaviour of a steadily de- 
creasing count rate over the course of the exposure is not 
in doubt. For the first 400 s l/-band exposure a higher res- 
olution light curve was created by splitting it into 4 s bins, 
and the second into 20 s bins, consistent with the majority 
of the exposures during that time. 

In order to get the best measurement of the optical 
temporal decay we created a single UV/optical light curve 
from all the UVOT filters. To do this the light curve in 
each filter was individually fitted to find the corresponding 
normalisation, and this was then used to renormalised each 
light curve to the V-band. The combined UV and optical 
light curve, from T-l-85 s up to T-l-2.5 x lO" s, is best fit 
with a power law decay with temporal index ao-pt = 1.64 ± 
0.01. The combined light curve shows no apparent colour 
evolution, illustrated in Fig. [2] 

The combined UVOT light curve is shown next to the 
0.3-10 keV and 15-350 keV light curve in Fig. [21 and a close 
up of the UVOT light curve from T-l-1000 s to T-|-30,000 s is 
shown in the inset. The log-log scale in this figure shows the 
decaying light curve in the 15-350 keV energy band from 
80-500 s much better than the linear scale in Fig [1] The 
magnitude along the left axis corresponds to the V^-band, 
and the count rate along the right axis applies to the X-ray 
and 7-ray light curves. 
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2.4 Observations By Other Facilities 

The prompt emission from GRB 061007 also triggered 
Konus-Wind and Suzaku, and the afterglow was monitored 
by ground based telescopes, which carried out both spec- 
troscopic and photometri c observations. GRB 061007 was 
measured to be at z=1.26 (|Osip et al.ll2006l : lJakobsson et al.l 
l2006l ). for which the rest-frame isotropic energy release and 
maximum luminosity in the 20 keV-10 M eV is Ejao ~ 1.0 x 
lOf^erg, and U^o ~ 1.8 x lO'^^ erg s"^ (|GoIenetskn et all 
I2OO6I ). The optical afterglow was also detected by ROTSE at 
26.4 s after the BAT t rigger at a magnitude of 13 .6 mag in 
an unfiltered exposure llRvkofffc Ruiopakarnl2"00^ . and the 
Faulkes Telescope South measured an fl-band magnitude o f 
R = 10. 15 ±0.34 142 s after the trigger (|Bersier et aLlbOOd ). 
consistent with the UVOT observations and indicative of a 
significant brightening during the first 90 s. The Swift obser- 
vations from T-l-80 s show that by this time the afterglow 
had entered a constant power law decay (see Fig. Ra- 
dio observations were also performed 1-1.24 days after the 
prompt emission, which resulted in no detection of a radio 
source, and a flux measure ment of 1 ± 36 fiJy at 19 GHz 
Ivan der Horst fc RolllioO^ . 

2.5 Multi- Wavelength Light Curve Analysis 

At ~ T + 100 s, after the last flare in the BAT hght curve 
(Fig. [1]), the BAT light curve appears to decay smoothly 
and monotonically out to ~T-(-500 s, after which the flux 
becomes indistinguishable from the background in the BAT 
energy range. A power law fit to the BAT light curve from 
T-l-lOO s to T-l-500 s provides an acceptable fit (x^ = 32 for 
29 dof) with a best-fit decay index of asAT ~ 1.61 ± 0.14, 
which is consistent with the decay observed in both the 
UV/optical and X-ray light curves. Furthermore, a power 
law spectral fit to this last part of the BAT light curve, 
beyond ~T-|-100 s, gives a best fit energy spectral index 
I3bat{> T + 100 s) = 0.80 ± 0.09, in fairly good agreement 
with the spectral indices determined for the X-ray and op- 
tical afterglows. This supports a scenario in which the 7- 
emission at this epoch is generated by the same radiation 
mechanism as the X-ray and UV/optical afterglow. 

2.6 The Spectral Energy Distribution 

Using the full range of data available with Swift an SED of 
the afterglow was produced using optical. X-ray and 7-ray 
data, which were normalised to an instantaneous epoch at 
600 s after the BAT trigger. At this time there are data avail- 
able from all three instruments, and good colour information 
from UVOT. For each of the UVOT lenticular filters the tool 
UVOt2pha (vl.l) was used to produce spectral files compat- 
ible with XSPEC (vl2.2.1) from the exposures taken closest 
to T+600 s. The count rate in each band was then nor- 
malised to the estimated value at T-f 600 s determined from 
the best-fit models to the light curves. The XRT spectrum 
was produced from WT mode data in the 0.5-10 keV energy 
range, omitting the first 30 s of data, which suffered from 
pile-up. This was scaled to the corresponding count rate at 
T-l-600 s. For the BAT spectrum, data in the 15-150 keV 
energy range were taken from the time interval T-l-lOO s - 
T+500 s, during which the BAT light curve decayed at a 
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Figure 3. The combined 7-ray, X-ray and UV/optical spectral energy distribution of GRB 061007 at T-l-600 s (see text) with best-fit 
models for each corresponding dust extinction curve shown; SMC (dashed), LMC (solid) and Galactic (dotted). Open circles are data 
points at wavelength A < 1215 A in the rest frame, and therefore not used in the spectral fitting. The rest-frame positions of the Galactic 
absorption feature at 2175 A the beginning of the Lyman-Q forest at 1215 A and the Lyman-a break at 912.5 A are indicated. 



rate consistent with X-ray and UV/optical afterglows and 
had a similar spectrum. This spectrum was normalised to 
T-l-600 s in the same way as the XRT spectrum. The SED 
is shown in Fig. [S] 

The lack of any break in the X-ray and UV/optical 
hght curves up to 10^ s places the cooling frequency ei- 
ther above the X-ray band or below the optical at T+600 s. 
Both energy bands will, therefore, have the same spectral 
index, and the SED was, thus, fit with a power law emis- 
sion spectrum, and two dust and gas components were in- 
cluded to model the Galactic and host galaxy photoelectric 
absorption and dust extinction. The column density and 
reddening in the first absorption system were fixed at the 
Galactic values (see sections 12.21 and 12. 3p . The second pho- 
toelectric absorption system was set to the redshift of the 
GRB and the neutral hydrogen column density in the host 
galaxy was determined assuming solar abundances. The de- 
pendence of dust extinction on wavelength in the GRB host 
galaxy was modelled on three extinction laws taken from 
observations of the Milky Way (MW), the Large Magellanic 
Cloud (LMC) and the Small Magellanic Cloud (SMC). The 
greatest differences observed in these extinction laws are in 
the amount of far UV (FUV) extinction, which is greatest 



in the SMC and least in the MW, and in the strength of the 
2175 A absorption feature, which is most prominent in the 
MW and negligible in the SMC . To model these extinction 
laws the parameterisations from lPeil (|l992l ) were used, where 
Rv = Av/E(B -V) = 3.08, 2.93, and 3.16 for the Galactic, 
SMC and LMC extinction laws respectively. 

An important difference between the Magellanic Clouds 
and the Milky Way is the average metallicity of these galax- 
ies, which is larger in the case of the Milky Way. The 
value that we determine in our spectral analysis for Nh is 
an equivalent neutral hydrogen column density that results 
from the amount of soft X-ray absorption in the spectrum, 
which will, therefore, depend on the abundances assumed in 
the spectral fit. The SMC and LMC have a metallicity that 
is ~ 1/8 and ~ 1/3 that of the Milky Way (|Peilll992l ). and 
this will result in an underestimation of the neutral hydro- 
gen column density by a factor of ~ 8 and ~ 3, respectively, 
if solar abundances are assumed. This is taken into account 
when we discuss the results from our spectral analysis in 
section 13.41 

At a redshift of z — 1.26 the beginning of the Lyman- 
a forest is redshifted to an observer-frame wavelength of 
~ 2580 A, which falls squarely on the UVWl fiher, the 
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Table 1. Power law spectral fit results to the spectral energy distribution of GRB 061007 



Spectral Model 


Host galaxy 
Nh^ (lO^i cm-2) 


Rest-Frame 

visual extinction, Ay 




(dof) 


Null- hypothesis 
probability 


SMC 
LMC 
MW 


O.DD_Q 2g 
r, or,+0.33 


0.39 ±0.01 
0.66 ±0.02 
1.13 ±0.10 


0.90 ±0.005 
0.98 ± 0.007 

-I -,0+0.04 
-^•-^'-'-0.03 


398 (308) 
347 (308) 
1042 (308) 


4 X 10-* 

0.063 

0.000 



Equivalent hydrogen column density assuming solar abundances at the GRB host galaxy. 



reddest of the UV filters. The unknown equivalent width 
of Lyman-a absorption from the host galaxy of the GRB 
and the variation of the Ly-a forest along different lines 
of sight makes it difficult to model this additional source of 
absorption in the SED. We, therefore, exclude the three data 
points affected by Lyman-a from our spectral analysis. The 
results from our spectral analysis are summarised in Table [T] 
The fitted value of the spectral index depends, for the most 
part, on the amount of optical/UV extinction, which in turn 
is affected by the dust extinction model used. Similarly, the 
amount of soft X-ray absorption fitted to the model will 
affect the best-fit spectral index. 

In all three cases an absorption and extinction system 
is measured at the GRB host galaxy, which is smallest in an 
SMC model and largest in a MW model. The lack of any 
evolution in Nh or (3 in the X-ray data (see section [23J, in- 
dicates that Ep has already moved through the X-ray band, 
and, therefore, the neutral hydrogen column density mea- 
sured is indicative of the amount of X-ray absorption at the 
GRB host galaxy rather than an artifact of evolution in the 
spectral index. The optical, X-ray and 7-ray spectrum is best 
fit by a model that assumes the GRB host galaxy to have an 
LMC extinction law for the GRB host galaxy, which gives 
St = 347 for 308 degrees of freedom (dof) with a null- 
hypothesis probability of p = 0.063. Models with an SMC 
and MW extinction law at the host galaxy were rejected 
with 99.96% confidence and 100% confidence, respectively. 



3 DISCUSSION 

GRB 061007 is one of the brightest bursts to be ob- 
served by Swift, with an isotropic energy rele ased in 7- 
rays that is second only to GRB 050904 (jTagliaferri et all 
I2OO5I ) and an optical flux compa rable to that observed in 
GRB 990123 (lAkerlof et al.lll999l ). The similarity of the af- 
terglow of GRB 061007 in the flux decay of 7-ray, X-ray and 
UV/optical bands, and the lack of features in X-ray and op- 
tical light curves distinguishes this burst from others. 

Optical afterglows typically decay as power laws, al- 
though usually with decay indices ~ 1.1 (Zeh, Klose & 
Kann 2006), and in the X-ray, the GRB hght curve is 
usually characterised by several discrete segments, with 
a canonical shape that is made up of three, and some- 
times four, power law segments with di fferent decay indices 
l|Nousek et all I2OO6I : IZhang et aLll2006l) . The temporal be- 
haviour of GRB 061007 is, therefore, not typical of other 
Swift GRBs, of which only ~ 7% observed promptly by the 
XRT (within the first few minutes of the BAT trigger) show 
a constant decay with no additional features such as flaring 
l|WiUingale et al.ll2006l ) (e.g. GRB 051117B, GRB 060403). 



When only considering long GRBs, the fraction of smooth 
decaying, featureless X-ray afterglows is even more infre- 
quent, amounting to only ~ 4% of the long GRB popula- 
tion. 

According to the fireball model, the temporal and spec- 
tral indices of the afterglow depend on the properties of both 
the GRB ejecta and the circumburst environment, which 
determine the location of the characteristic synchrotron fre- 
quency, Vm and the cooling frequency, t^c, with respect to 
the observing band l|Zhang et al.|[2007al ). At any one time 
the location of these frequencies are dependent on the energy 
distribution of electrons in the circumburst environment de- 
fined by the energy index p, the density and density profile 
of the circumburst medium, the amount of kinetic energy 
in the outflow, and on the fraction of energy in electrons 
and in the magnetic field, te and es, respectively, and 
e_B are typically assumed to remain constant with time, and 
therefore the evolution of Vm and Vc, and thus the afterglow 
decay and spectral indices, depend on the density and den- 
sity profile of the circumburst medium and on the amount 
of kinetic energy released by the GRB. 

The lack of a break in the light curves provide 
constraints on the spectral regime that the X-ray and 
UV/optical afterglow are in, and on the frequency of fm 
and Vc, neither of which could have crossed the observing 
window from T+80 s to T+10^ s. From the white filter we 
determine that around half the counts detected in the white 
band observations correspond to the V and B filter. The 
lack of a break in the first white filter exposure, therefore, 
indicates that at T+80 s, Vm was no bluer than the B-band 
filter. This corresponds to an upper limit on the character- 
istic synchrotron frequency Um < 7.5 x lO'^* Hz. From the 
UV/optical data we also determine that at T+80 s the flux 
density in the V-band, corresponding to v — 5.5 x 10^'* Hz, 
is — 730 mjy, after correcting for absorption at the host 
galaxy. 

With these constraints in place we investigate the var- 
ious conditions that can explain the temporal and spectral 
properties of GRB 061007, all of which can be categorised 
into two parent models that differ fundamentally in the ef- 
fect of coUimation on the observations. In the first model we 
assume no jet break to have occurred before T+2.65 x 10"" s, 
based on the constraints we obtained in section 12.21 and 
the observed emission is, therefore, indistinguishable from 
isotropic emission. Secondly, we investigate a model in which 
the jet break occurs before our first X-ray and UV/optical 
observations at T+80 s such that the assumption of isotropic 
emission no longer applies. 
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3.1 Spherical Expansion Model 

The simplest model to consider is one in which the observed 
afterglow emission follows the flux decay and spectral evolu- 
tion as would be expected from synchrotron radiation emit- 
ted from an isotropic outflow. As already mentioned, one 
of the important factors that determines the evolution of 
the defining synchrotron frequencies is the density profile 
in the surrounding circumburst material, and we therefore 
consider both a wind-like circumburst environment, such as 
would be expected in the vicinity of a massive star, and a 
constant density profile, such as in the ISM. For simplicity, 
we assume the density profile of the circumburst medium to 
have the form p oc r"" , where s = and s — 2 corresponds 
to an ISM and wind-like profile, respectively. 

The relation between a and /? for synchrotron emis- 
sion at different intervals in the synchrotron spectrum are 
defined by a set of closure relations, both for a constant 
den sity medium (ISari et al.i il998.) and a wind density pro- 
file (IChevalier &: LilboOOl ). Assuming spherically equivalent 
emission, the standard closure relations do not fully satisfy 
the observed spectral and temporal indices of GRB 061007 
for either an ISM or a wind-like circumburst medium, where 
Qo,x = 1.65 ± 0.01 and f3o,x = 0.99 ± 0.02. This could be 
the result of additional components unaccounted for in the 
closure relations, or an oversimplification of the assumptions 
made when outlying the closure relations e.g. the density 
profile may not be distributed as r'". 

For both a wind and an ISM circumburst density pro- 
file our observations are in closest agreement with a model 
in which the observed synchrotron emission is in the slow 
cooling phase {vm < Vc), where the UV/optical and X-ray 
bands are in the same slow-cooling regime {vrn < v < Vc). 
The decay index is Aa 0.17 steeper than expected for a 
constant circumburst density profile, and requires a fairly 
soft electron distribution index p ~ 3. In the case of a wind 
circumburst environment the decay index is shallower than 
predicted by this model, by a difference Aa ~ 0.33. It is pos- 
sible that the density profile in the environment surrounding 
GRB 061007 lies somewhere between the two scenarios con- 
sidered here, where < s < 2. Alternatively, the shallowness 
of the decay index for a wind case could be explained if en- 
ergy injection is involved. This would require a late-time 
source of energy into the forward shock that would slow 
down the cooling rate of the excited electrons and, thus, 
slow down the afterglow decay. Either way, the lack of a 
cooling break for the duration of Swift observations out to 
10^ s requires consideration, and the conditions necessary to 
maintain fc above or below the observing frequency range 
need to be investigated. 

The evolution of i^c with time is dependent on the den- 
sity profile of the circumburst medium, whereby for s > | 
Vc will increase w ith time, whereas for s < | it will decrease 
(|Panaitescdl200(j ) . The conditions needed for the cooling fre- 
quency to not have crossed the observing window, therefore, 
differ for an ISM and wind circumburst density profile. For 
an ISM medium the cooling frequency has to be maintained 
above the X-ray band for up to T-|-10^ s, whereas in a wind 
scenario the cooling frequency needs to be above the X-ray 
band at the start of the X-ray observations, at T-l-80 s. Ei- 
ther way the cooling frequency must lie above the peak X-ray 
frequency out to at least T-|-10^ s for a cooling break not 



to cover in the X-ray observations (see section I2.2|l . This 
corresponds to the requirement that i^c > 6 x 10^^ Hz at 
T-|-10^ s. In the next section we investigate the consistency 
of a spherical emission model with the data for both an ISM 
and wind medium where Vrn < < Vc- 

3.1.1 ISM Case 

For an ISM-like circumburst environment, the temporal and 
spectral fiux decay indices are related as a = 3/3/2. For a 
value of Px = 0.99 ± 0.02, we expect a = 1.49 ± 0.03, which 
is shallower than the observed decay rate of a = 1.65 ±0.01. 
However, this deviation from the relation is small, and could 
be accommodated if the density profile is slightly different 
from the uniform density we have assumed. In this regime 
P — {p ~ l)/2, where p is the electron energy index. Thus, 
for l3x = 0.99±0.02 the electron energy index is in the range 
2.94 < p < 3.02. We therefore take p = 3.0. These relations 
hold as long as the characteristic synchrotron frequency, Vm , 
lies below the optical band, i.e. Vm < 7.5 x 10^* Hz, at 
T-fSO s, as discussed previously. 

For a constant density medium in which synchrotron 
radiation is the dominating emission mechanism, the char- 
acteristic synchrotron frequency, the cooling frequency and 
the peak density flux are given by the following equations 

o o ^, inl5 /'P ~ , n1/2 1/2 2 ^1/2 ,-3/2 -rj 

i^,„ = 3.3x10 \^—jj{l + z)'ej^ ee£fe;52*d Hz (1) 
= 6.3xlO''(l+^)-'/'(l+F)-'e~'/'£;-5f n-'t;'/' Hz(2) 

F.,max = 16(1 + z)D28ey^Ek,52n^^ mJy (3) 

jZhang et al.ll2007al '). where es is the fraction of energy in 
the circumburst magnetic field, ee is the fraction of energy 
in electrons, -Efc,52 is the kinetic energy in the GRB in units 
of 10^^, td is the time since the prompt emission in units 
of days, n is the density of the surrounding circumburst 
medium, Y is the Compton parameter, and D2S is the lu- 
minosity distance in units of 10^* cm. Typical values for the 
fraction of energy in the cicumburst magnetic field and in 
the electrons are es = 0.01 and ee = 0.1, respectively. 

For a spectral index /3 = 1, vF,^ is constant for i/m < 
u < Uc. At T+80 s = 730 mJy at u ^ 5.5 x lO^"* Hz, 
and therefore i^mfi^.max ~ 730 x (5.5 x lO^**). By using this 
together with Eq. [T] and |3] we obtain: 

3.4x10"'^ 

" = ,2 ,4^3 (4) 

In Eq. [2]we apply the condition that at 10^ s after the 
prompt emission Uc must still be greater than 6 x 10^^ Hz, 
and substituting in Eq.|4] gives 

(i + y)-^4/'4<52>o.52 (5) 

Considering now the characteristic synchrotron fre- 
quency, we know that Um< 7.5 x 10^* Hz, and this reduces 
Eq.[T]to 

^T^lEl%< 1.7x10-', (6) 

This, and the inequality given in Eq. [S] provide an upper 
limit to the value of e_g in terms of -E^ 52, such that 

es < 3.1 x 10-1^(1 + y)-*£^,52 (7) 
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The above equation bounds us to extreme values of Ek^52 
and eg- For Ek,z2 = 10^ es < 3.1 x 10""'(1 + Y)-^, which 
is an unrealisticaUy low value of es. Even a value as small 
as es = 10~® requires a huge amount of kinetic energy to 
meet the condition that Ek,52 > 1-5 x 10*(l + y)'*/^ erg. The 
amount of kinetic energy assumed is an isotropic-equivalent 
value which would, thus, be smaller if it were corrected for 



coUimation. A 3a lower limit of ti 



2.65 X 10'' s can be 



imposed as the earliest time that a jet break could have 
occurred, tj is the time of the jet break measured from the 
onset of the afterglow, which we take to be equal to T. A 
kinetic energy Eh,52 = 1.5 x 10^(1 + F)"''^ erg and es = 
lO"** would require = 0.012(1 + F)"^/^ and n = 49(1 + 
Y)~^^^ cm~"^ for the conditions given above to be satisfied. 
From this we estimate the jet opening angle using 



1/8 



+ ZJ \E ,,,2J ' 

(jPanaitescu fc Kumar! [2OO0I ). where tj,d is in units of days. 
The equation for the IC parameter is as follows: 

l + (l+4»7ir?2ee/eB)'^^]/2 (9) 

(2-p)/2l 



Y 



Zhang et al.l l2007al ). where t]i — min[l, {ydv-n., j 
ban fc EsinI [200li r and <C 1 accounts for the Klein- 



Nishina correction. For the physical parameters assumed, 
the IC parameter, K, is near negligible (^^(1 + F) < 
2.2 X 10""), and can therefore by approximated to y = 0. 
In this case, a jet opening angle for GRB 061007 would 
need to be Qj > 6.3° to produce a jet break later than 
T+2.65 X 10* s. This would provide a beam-corrected kinetic 
energy of E^"^'' > 9 x 10*^ erg, which is half a solar mass 
and three orders of magnitu de greater than usual bursts 
(IFreedman fc Waxman|[200ll ). The energetics required by 
this model are unreasonably large and, therefore, an unlikely 
solution to explain the properties of GRB 061007. 



3.1.2 Wind Case 

In the slow cooling regime where Vrn < v < Vc, the spectral 
and temporal indices for a wind-like circumburst environ- 
ment are related such that a = (3/3 + l)/2, which can be 
rejected by the observations at more than 10 a confidence. 
In order for this model to be compatible with the data an 
additional source of energy is required to slow down the 
decay rate of the afterglow. A continual source of energy in- 
jection could come from either the central engine itself, or 
late time refreshed shocks that slow down the decay of the 
afterglow. The shallow decay phase observed in many XRT 
light curves is interpreted to be the result of energy injec- 
tion, and although this phase typically lasts for 10* s, there 
are also examples of GRBs with evidence of energy injection 
for prolonged periods of time, and with very smooth decay 
slopes (e.g. Romano et al. 2006; Huang, Cheng fc Cao 2006). 
Taking this into account, an explanation whereby continual 
energy injection maintains the afterglow and slows down its 
decay over a long period of time is not unreasonable. For the 
purpose of this analysis the source of this additional energy 
is not considered and the energy injection is simply quanti- 
fied by the expression L oc t~'', whe re L is luminosity in units 
of erg s~^ (|Rees fc Meszaroj|l998l ). A wind circumburst en- 
vironment may then be compatible with our observations 
ii < V < Vc and q = 0.66 ± 0.02 and p = 2.98 ± 0.02. 



Using these best fit parameters and the known constraints 
provided by the data it is possible to determine the param- 
eter space of the GRB kinetic energy and the circumburst 
medium microphysical parameters that could reproduce the 
observations. 

For a wind environment the characteristic synchrotron 
frequency, i/m, and flux at this frequency are given by 



4.0 X 10 



1/2 ,-3/2 TT 



{p-0.69){l + zy^^ 



(10) 



77{p+0.12){l+zf^D^ie'J^E]ll^A,t-'^^m.Mn) 



where A, — A/ (5 x 10^^) g cm ""^ is t he density scaling such 
that p = Ar-^ (| Chevalier fc Lill2000l '). As before the flux at 
u = 5.5 X 10" Hz is F^ = 730 mjy at T-l-80 s, and since 



A, = 



I^m-Fi/,max, WC haVC 

9.16 X lO"'^ 



eseeEK, 52 



(12) 



The dependencies on the cooling frequency also vary from 
Eq. [2]for a wind medium, where it is now defined by 



V, = 4.4 X 1010(3.45 -p)eO'«P(l -I- ^)-3/2 



Yr'^B'''E^ilt2A7X^' Hz, 



(13) 



which has been modified from IChevalier fc H (l2000fl to in- 
clude the Compton parameter, Y . A new expression, there- 
fore, needs to be defined regarding the conditions on the 
cooling frequency. Since i/c increases with time, the only 
useful constraint is a lower limit on i^c at early times. This 
constraint can be determined through spectral modelling of 
the afterglow SED at an epoch of T-l-300 s, which is the 
weighted mid-time of the time interval over which the BAT 
spectrum shown in Fig. [3] was accumulated. This SED was 
fit with a broken power law with a spectral break above 
the BAT energy range, and the change in spectral index at 
the break was fixed to A/3 = 0.5 to correspond to a cooling 
break. Taking the from this fit into account, the Vc upper 
limit was then determined by shifting the spectral break to 
lower energies until the goodness of fit deviated by Ax^ ~ 9 
(3(7 confidence for one interesting parameter). This gave an 
upper limit on the cooling frequency of Uc > 1.2 x 10^^ Hz 
at T+300 s. Applying this, Eq. [T3] combined with Eq. [TD 
reduces to 



(1 + Y)-\]^\tE''Jl2 > 0.0073 



(14) 



As for an ISM circumburst environment, Um < 7.5 x 10^*, 
which if applied to Eq. [10] gives 



(15) 



1/2 2 j-,1/2 ^ e DC „ in-6 
eg Ee-Ej^ 52 < 6.26 X 10 , 

and combining this expression with Eq. [TJ] gives 

eg < 2.88 X 10''^(1 + y)"*£K,52 (16) 

For e_B = 10~® the kinetic energy in the forward shock 
needs to be Ek,52 > 3 x 10^(1 -|- Y)^ erg, which when 
taking into account the possible jet coUimation reduces to 
El^"' > 1.8 X 10'^^(1 -I- Y)* erg. However, it is also neces- 
sary to include the energy injected into the fireball, which 
increases the total energy in the forward shock with time 
as Ek oc t^~'^ . For q = 0.66 the total energy in the for- 
ward shock at T-l-10* s increases by an order of magnitude 
to S^"'"'' > 1.8 X 10*''(1 + Yf erg, which is similar to the 



10 p. Schady, et al. 



amount of kinetic energy required for an ISM-like circum- 
burst environment. A wind environment, therefore, faces the 
same energy budget problem as present for an ISM environ- 
ment 

In the spherical emission model discussed above the pri- 
mary assumptions made in both the ISM and wind scenarios 
are that both the X-ray and UV/optical afterglow observed 
are produced by the same radiation mechanism, and that 
the decay rate is equivalent to spherical emission. The dif- 
ficulty in accommodating the observations with a standard 
fireball model indicate that in the case of GRB 061007 either 
synchrotron emission is not the dominant radiation mech- 
anism throughout the observing energy range, or that the 
GRB outfiow is coUimated whereby the observed emission 
is no longer equivalent to isotropic emission. With this in 
mind we investigate how additional radiation mechanisms 
and various collimation effects would alter the observed af- 
terglow from that produced by synchrotron emission from 
an isotropic outfiow. 



3.2 Highly Collimated Models 

The large amount of energy required to power a model with 
isotropic emission leads us to investigate a scenario in which 
the emission from the GRB is highly collimated, and we once 
again assume the circumburst environment to be of constant 
density. A very narrow jet opening angle could reduce the 
energy required to explain our observations by several orders 
of magnitude, where tj must be prior to our first afterglow 
observations in order to account for the lack of a break in 
the X-ray and UV/optical light curves at later times. If we 
assume the afterglow begins no earlier than the time of trig- 
ger, this corresponds to a time tj < 80 s. Within the premise 
that the outflow is highly collimated we consider two cases in 
which the cooling frequency is either above or below the ob- 
serving frequencies. In these two regimes p and /3 are related 
hy P — {p ~ l)/2 and /3 = p/2, respectively, corresponding 
to p ~ 3.0 or p ~ 2.0. Aft er the jet break the flux decays as 
Ft oc t-P l|Sari et al.lll999l ). where p is either 2 or 3, which is 
steeper than the observed decay rate of a ~ 1.65. However, 
this does not account for continual energy injection, which 
would reduce the afterglow decay rate. 



3.2.1 Vm< V < Vc 

For a laterally expanding jet with constant density circum- 
burst medium and v < Vc, the relation between the 
temporal decay of the GRB afterglow, a, and the spectral 
index, /3, post jet-break is given by 



oc t 



3/8 



a = (l + 2/3)--(l-q)(/? + 2). 



(17) 



ijPanaitescu et al] l2006l ) where q is t he energy injection 
param eter such that L = Lo{t/to)~'' (jZhang fc Meszaroj 
I2OOII ). and since we take the time of our jet break to be 
at tj = 80 s, to = tj, L = Lo{t/tj)~'' . For our observed 
temporal and spectral indices, q = 0.3 and p = 3, as be- 
fore. Prior to the jet break Eq. IH3I still apply, but in or- 
der take into account the effect of energy injection, the ki- 
netic energy requires an additional factor of (t/fj)'^"'''. The 
time of jet break and the jet opening angle are related as 



To keep 9j as large as our observations will al- 
low we take tj = 80 s. At T-(-80 s the additional factor of 



\(l-9) 



is unity and there is no change in Eq. [T]|3l Thus 



at T-l-80 s, Eq. gland Eq. |S]stand. For Eq.[2]the constraint 
that Uc > 6 X 10^'' at T-|-10^ s means that the evolution of 
Uc post jet-break has to be taken into account. If Uc^ is the 
frequency of Vc at the time of the jet break, T-(-80 s. 



: 4.05 X lO^^l + Y) 



-2 1/2 4 p5/2 



(18) 



However, we are concerned with where Vc lies at T-l-10 s. 
After a jet break oc 7~*t~^, where 7 oc ^"(^+5)/^ is the 
Lorent z factor of the electron s in the circumburst environ- 
ment (jPanaitescu et ai1l2006l ). The cooling frequency after 
T-l-80 s is given by 

where 7j is the electron Lorentz factor at T+80 s. Using the 
condition that Uc > 6 x 10^^ Hz at T+IO'"' s, for q = 0.3 
Eq. [19] reduces to 



{l + Y)-\y'etEl%>0.41 



(20) 



In the same way that we determined the inequality given 
in Eq. [3 we combine the above equation with Eq. [S] and 
determine that 



es < 5.0 X 10~''(l-Hy)-*£^,52 



(21) 



For Bfc.sa = 5 x 10^ erg, eg < 6.2 x 10"''(1 + Y)-\ 
and if we take this upper limit as the value of es, te ~ 
3.1 X 10~^(1 + Vy^^. Although the assumed kinetic en- 
ergy is large, this corresponds to the isotropic-equivalent 
value and the beam-corrected kinetic energy is a few or- 
ders of magnitude smaller than this. For p = 3, the in- 
verse Compton parameter Y is proportional to [r;(ee/e_B)]^''^, 
where ri=mm[l, (um/iyc)^^^], and for these values of te and 
es, F < 0.3. For values, of £^,52 = 5x lO", es = 6.2 x 10"^ 
and Ee = 3.1 X 10~^ the particle density in the circumburst 
environment is n = 0.08 cm~'^, which, using Eq. [8] gives 
an opening angle of 6j = 2.0 x 10^"^ rad. This corresponds 
to 0j = 0.11°, which gives a beam-corrected kinetic energy 
^corr ^ ^q5i ^^.g energy = 2x 10*® erg. This 

would imply an efficiency of 77 ~ 0.002, which is two orders 
of mag nitude smaller than is typically inferred (jZhang et al.l 
l2007al '). The jet opening angle that we estimate assumes 
the hydrodynamics as described by IPa naitescu &: Kumar! 
(2000), and variations on th e hydrodynamics, such as that 
discussed in lSari et al.l (|l999l ). result in a different coefficient 
in Eq. |S] However, these differences are of order unity, and 
do not affect the validity of the model discussed. 

This model is, therefore, able to satisfy our observations 
without the requirement of an excessively large kinetic en- 
ergy. However, the opening angle i s a factor of at lea st 30 
narrower than previously observed (iBloom et al.|[2003l ) and 
the inferred 7-ray energy is extremely low for a long GRB, 
prompting us to consider other options. 



3.2.2 v>vc 

Keeping in place the condition that tj < 80 s, we investigate 
a model in which the X-ray and UV/optical bands are above 
and Vc. In this model our observational constraints on 
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the values of Um, Vc and -Fi^.max only apply at T+80 s. Since 
we assume a jet break to occur at tj = 80 s, the factor of 
(t/tj)'^"''-' in the kinetic energy is unity, and q is therefore 
not of importance in our calculations. For Uc above the ob- 
serving band, (5 = p/2, and therefore p = 2, in which case 
the dependence on the minimum Lorentz factor of the shock 
accelerated electrons 7m changes from 



_ (p-2) mp 



(p - 1) mo 



to 



7m 



In 



(-) 



6 

nie 



(22) 



(23) 



where Sm and Em are the minimum and maximum energy 
of the shock accelerated electro ns, and F is the bulk Lorentz 
factor in the shocked medium l|Pe'er fc Zhand [20061 ). Vm oc 
7m , and therefore the expression for Vm given in Eq. [T] has 
to be modified such that [{p — 2)/{p — 1)]^ is replaced by 
[ln(£m/£Af )]^. The value of ln(eM /cm) is not well determined 
due to the unknown magnetic field in the upstream of the 
GRB shock, although it generally ranges between 5 and 10 
l|Li fc WaxmanlbOOd ). We therefore take ln(£Af/£m) = 7, 
such that Eq. [T] becomes 

i.m = 6.7 X 10" Hz(l + z)'^^e'J^e^,El%tf^^ (24) 

We assume the same jet break time as before {tj = 80 s), 
whereby Eq. [5] and [53] reduce to 



and 

^'4eI^^. < 2.1 X 10- 

In the slow cooling regime, where < J^c, 

/ \ -(p-2)/2 / N -p/2 



(25) 
(26) 

(27) 



By taking the measured value of Fi, — 730 at = 5.5 x 
10^^ Hz, and using Eq. (2] [3] and 1241 the expression for F,^ 
reduces to 



ee£;fc,52(l + F)"' = 3.65 



(28) 



From these equalities we determine that a kinetic energy 
of Ek,s2 = 10^ erg would require ee = 3.7 x 10"^(1 + Y), 
€B < 0.25(1 + y)-* and n > 0.05(1 + Y)^. For a GRB local 
environment with a particle density of n = 10'^ cm""^, which 
is a t ypical density for a molecular cloud (IReichart fc Pric"3 
|2002| ). the jet opening angle would be 6j = 0.8°, a factor 
of a few smaller th an that determined f or some other GRBs 
(e.g. GRB 980519; iNicastro et"al]|l999l ). This gives a beam- 
corrected kinetic energy of E^°^2 = 9.7 x 10®° erg. 

This model provides an alternative scenario in which a 
highly coUimated jet is the cause for the large fiux observed 
in GRB 061007, but where a larger jet opening angle and 
more standard beam-corrected 7-ray can satisfy our obser- 
vations, with respect to the previous model. 



3.3 Synchrotron Self-Compton Emission 

For completeness we look into a scenario in which syn- 
chrotron self-Compton (SSC) emission contributes to the 



observed emission. In the case where the X-ray emission 
is dominated by SSC radiation, and the UV/optical emis- 
sion is produced by synchrotron radiation, both observed 
bands would need to be in the fast-cooling regime, with 
I'opt > > Vm and ux > I'i'^ > v^^' , in order for the 
decay slope and spectral index in the X-ray and UV/optical 
emission to be consistent. In this case the spectral index is 
given by f3 = p/2, which would give a value of p = 2.0. 
However, the decay of the optical afterglow is expected to 
be —(2 — 3p)/4, which is inconsistent with our observations 
of a = 1.65 ±0.01. 



3.4 The Circumburst Medium 

From our spectral analysis we determine the dust extinction 
in the circumburst medium to be most consistent with a 
wavelength dependence similar to that of the LMC. This is 
consistent with previous studies of GRB environments (e.g. 
Fruchter et al. 2006, Kann, Klose & Zeh 2006; Schady et al. 
2006), which also find GRB hosts to be irregular galaxies, 
and is also in agreement with the coUapsar model, which re- 
quires a sub-solar metallicity progenitor star, with an upper 
limit of Zq < 0.3 (e.g. Hirschi, Meynet & Maeder 2005; 
Woosley fc Heger 2006). 

The best-fit parameters for the X-ray column density 
and dust extinction local to the GRB are Nh ~ (5.3 ± 
0.3) X 10^^ cm"^ and Ay = 0.66 ± 0.02, and this gives 
a gas-to-dust ratio of Nh/Av = (8.0 ± 0.5) x 10^^ cm"^. 
However, Ay was determined assuming an LMC dust ex- 
tinction law, whereas Nh was determined assuming solar 
abundances, which are unlikely to be correct for an irreg- 
ular, Magellanic-type galaxy. Correcting the X-ray column 
measurement to a metallicity of 1/3 solar, which is appro- 
priate for an ISM similar to that of the LMC, we find that 
the gas-to-dust ratio is a factor of more than three larger 
than that measured in the LMC. We note that an environ- 
ment's metallicity affects the amount of dust present as well 
as the equivalent Nh measured. Our assumptions on the 
metallicity, therefore, do not change the Ay / Nh ratio as we 
measure it, since to first order both Ay and equivalent Nh 
scale directly with metallicity. It is, therefore, valid to com- 
pare the dust-to-gas ratio in different environments derived 
this way, independent of their metallicity. Although not con- 
sistent with the LMC, the dust-to-gas ratio measured in the 
local environment of GRB 061007 is consistent with the av- 
erage gas-to-dust measured in the local environme nt of a 
sample of previous Swtft GRBs (|Schadv et al.ll2006l ). which 
was (Nn/Ay) = (6.7 ± 1.1) x 10^^ cm'^ 

The lack of evolution observed in the soft X-ray absorp- 
tion later than ~ 80 s after the prompt emission suggests 
that the burst must have completely ionised the gas in its 
immediate environment within this time scale. The absorb- 
ing medium that we observe must, therefore, be far enough 
from the burst that it was not completely ionised in the first 
few 100 s of seconds after the initial explosion. The probabil- 
ity of there being an absorption system in the line of sight 
from an unrelated, intervening galaxy that is of sufficient 
optical depth to contribute a measurable amount of soft X- 
ray absorption is ver y small (|0'Flahertv fc Jakobsenlll997l : 
ICampana et al.ll2006l ). leading us to conclude that the ab- 
sorbing material measured in the spectru m of GRB 061007 
is local to the burst. From the results in IPerna fc Lazzatil 
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l|2002l ) we estimate the absorbing medium to lie at a few 
tens of parsecs from the GRB, and the large value of Nh 
suggests that GRB 061007 was embedded in a dense, star- 
forming region, or giant molecular cloud, in common with 
many other GRBs (e.g. Dai & Lu 1999; Wang, Dai & Lu 
2000; Reichart & Price 2002). 



3.5 Comparing Between Models 

It is interesting to note that with Ej^iso = 10^"* erg in the 
10 keV-20 MeV energy range, and a rest-frame peak en- 
ergy of Ep = 902 keV, GRB 061007 satisfies th e Amati 
relation, given by Ep oc El^^f^^ (| Amati et al.ll2002l 'l. It also 
satisfies a correlation between the peak isotropic luminos- 
ity, Lp^iso, and spectral lag in the prompt emission, nag, 
whereby Lp^iso oc Tj"^'^. Both these correlations are sat- 
isfied by a large fr action of long GRBs (e.g. lAmatil [20061 : 
iGehrels et al] 120061 '). suggesting that GRB 061007 has the 
same class of progenitor as other long GRBs. 

The consistency observed in the spectral and tempo- 
ral afterglow behaviour in the X-ray and UV/optical bands, 
and in the late time 7-ray band, implies that the emis- 
sion observed comes from the same spectral segment. If we 
take the afterglow to be in the slow cooling regime and 
assume a jet break to occur at tj > 2.65 x 10'' s a large 
amount of kinetic energy, of at least E^°^'~ = 9 x 10^^ erg 
is required for a constant density circumburst medium, or 
j^corr y g Q ^ -1^052 ^ wiud-like circumburst medium. 

The implausibly large amount of energy required by both 
these scenario makes them unsatisfactory models to explain 
the properties of GRB 061007. 

We also investigate a second scenario in which a highly 
coUimated outflow produces a jet break before T-f 80 s. This 
requires a jet angle no larger than 0.8°, which satisfies the 
observational requirements with a lower energy budget and 
more reasonable set of parameters. The upper limit on the 
jet opening angle of 0.8° is a factor of 4 smaller than that 
previously determined in GRBs (Berger, Kulkarni & Frail 
2003), although still physically possible. Furthermore, the 
need to observe the jet break in order to measure 6j intro- 
duce selection effects, and the smaller and larger end of the 
current jet opening angle distribution will, thus, be under- 
estimated. For the very smallest of jet opening angles to be 
measured, very early time data of the afterglow is needed, 
preferentially covering several energy bands. 

For a jet break earlier than T-l-80 s our observations 
are valid for a spectral regime in which i/m < 1^ < Vc, and 
Vc < V, where in the latter a larger jet opening angle is 
obtained by assuming a density n = 10* cm~^, consistent 
with that observed in molecular clouds. The dependence of 
the circumstellar particle density on e_B, ee and Ek,52 in the 
former of these models does not allow for such high densities, 
for which we estimate n = 0.08 cm~^. Our spectral analysis 
on the afterglow resulted in a large column density, which 
supports a scenario in which the GRB is embedded within 
a dense molecular cloud or active star- forming region. This 
therefore favours a model in which Vc < v, where the number 
density is inferred to be larger. 

In the case of a uniform jet, where the energy and 
Lorentz factor are constant across the jet, the afterglow after 
a jet break should decay as t~^. In both the post jet-break 



models that we consider the predicted decay slope is steeper 
than the observed value of q = 1.65±0.01. A possible expla- 
nation is continual energy injection, which would maintain 
the energy in the afterglow and slow down the rate of de- 
cay. The requirement for steady energy injection over the 
entire light curve is of some concern, but seems plausible 
on the basis of many other examples of energy injection in 
Swift afterglow observations. Evidence of prolonged energy 
injection in the form of plateaus in X-ray light curves is 
observed in a large fraction of Swift GRBs, although typi- 
cally up to lO'^-lO* s after the prompt emission. In the case 
of GRB 061007 a steady energy injection is needed for at 
least an order of magnitude longer than this, during which 
no break is observed in the light curve. This is not to say 
that the central engine needs to be active for this length 
of time, but simply that it is active for the duration of the 
prompt emission and ejects shells with a large distribution 
in the Lorentz factor. A range in F can provide a source 
of energy to the afterglow for a long time after the cen- 
tral engine has stopped thr o ugh ongoing refreshed shocks 
(e.g. iRees fc MeszarosI Il998l : ISari fc MeszarosI |2000D . Nev- 
ertheless, an energy injection mechanism that remains ex- 
tremely constant for such a long time requires specific con- 
ditions to be in place, such as a central engine luminosity 
that varies smoothly with time in the case of a long-lasting 
engine, or a smooth distrib ution of Loren tz factors in the 
case of refreshed shocks ( Zhang et al.ll200Q) . 

Alternatively, the slow decay rate of the afterglow in the 
context of post jet-break could also be the result of a jet- 
edge effect in which there is no sideways expansion. In such 
a case the afterglow decay index would steepen by 3/4 when 
r = l/6j, and would conti nue to decay at this rate due to the 
lack of lateral expansion eszaros fc Reed 119991 ). In such 
a scenario there is no need for energy injection, although 
additional conditions would be required to prevent the jet 
from expanding further, such as a high magnetic field within 
the jet. 

A possibly more natural form of non-laterally expanding 
jet is a structured jet or outfiow with an energy distribution 
that is a decreasing function of angular distance from the jet 
axis, 60, and a narrow, high-energy core at the axis. A jet 
break is, then, observed when the edge of the uniform inner 
core becomes visible to the observer, and it is, therefore, the 
viewing angle rather than the jet half-opening angle that 
determines the time of the jet break. In order for a jet break 
to occur before T+80 s the viewing angle needs to be close 
to the jet axis, which is not that unreasonable given the 
brightness of GRB 061007. At the time of jet break the 
smoother distribution in energy, e, compared to a uniform 
jet results in a smaller change in the decay slope that will 
be shallower than the 3/4 increase in decay index observed 
in a uniform jet break. For a power law distribution of E, 
where E oc 9~'', the smaller the value of k the shallower 
the afterglow light curve will be after the jet-break. In the 
case of Vm < Vobs < I'c, where p = 3, a decay rate after a 
jet break with index a = 1.65 will occur in structured jet 
with an E distri but ion index k = 0.3, and for Pc. < i^ots we 
require fc = 1.4 (jPanaitcscu 200S). 

A consequence of a very early jet break time is the 
faster evolution of F into a non-relativistic phase, which 
will produce an additional break in the light curve as the 
decay rate slows down. For a uniform jet with no energy 
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injection, this brea k would occur at at a time t-nr = tjO- 
l|Livio fc Waxmanll200dl '). For tj = 80 s and Oj < 0.8°, the 
jet would, therefore, not enter the non-relativistic phase un- 
til more than 4 x 10^ s after the prompt emission, which is 
not in conflict with our observations, and in the presence 
of energy injection or for structured jets t„,. would only get 
larger. 

In the context of post-jet break emission, a decay slope 
of a = 1.65 can, therefore, be accounted for by a non- 
laterally expanding jet as well as continual energy injection 
through refreshed shocks. We, therefore, find an early jet 
break before T-l-80 s can satisfy the conditions imposed by 
our observations without requiring a large energy budget 
and with a reasonable set of physical parameters. 



4 CONCLUSIONS 

In this paper we have presented multi-wavelength spectral 
and temporal data taken with Swift of GRB 061007, an 
very bright GRB with unusually smooth panchromatic de- 
cay. Using these data we explored the circumstances that 
could produce such a bright burst, and so determined the 
conditions that distinguish GRB 061007 from other, more 
standard GRBs. On many levels we find GRB 061007 to be 
typical of other GRBs, suggesting that the same progenitor 
model for most long GRBs is also applicable to GRB 061007. 
The surrounding dust has properties most consistent with 
an LMC-like environment, and the circumburst material is 
likely to be dense, suggesting that it is embedded in a large 
molecular cloud or a region of active star formation. Fur- 
thermore, its consistency with both the Amati and timelag- 
luminosity relations, which apply to a large fraction of long 
GRBs, indicating that the source of energy in GRB 061007 
is no different than that of most other long GRBs. The enor- 
mous brightness of GRB 061007 is, therefore, unlikely to be 
due to a different type of progenitor, but the result of ex- 
treme values for a few defining characteristics. 

Two inherently different models can satisfy our obser- 
vations of GRB 061007, where the main difference lies in 
the time of the jet break; either before or after our obser- 
vations. Both require certain extreme conditions to account 
for the brightness of this burst, either in the form of a large 
energy budget or in the coUimation of ejecta into a very 
narrow jet opening angle. However, given the brightness of 
GRB 061007, it should not be a surprise that some of the 
defining properties are non-standard. The large amount of 
kinetic energy required in a late-time jet break model lead us 
to favour a model in which the GRB 061007 is produced by 
the concentration of energy into a highly narrow jet < 0.8°; 
a condition that can more easily be met than the production 
of large amounts of energy. 

The brightness of GRB 061007 across the electromag- 
netic spectrum makes it ideal to explore the conditions 
surrounding the production of GRBs, and its uniqueness 
provides further insight into the range in properties that 
GRBs have. The small opening angle that we determine for 
GRB 061007 has important implications on the true range 
in jet opening angles, and this could provide a clue as to 
the cause for the large fraction of Swift GRBs with no clear 
signs of jet breaks. 
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